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TECHNICAL NOTE NO. II 5 S p ^ 

FLUTTER AND OSCILLATING AIR-FORCE CALCtTLATIONS FOR AN 

AIRFOIL IN A TWO-DIMENSIONAL SUPERSONIC FLOW 

By I. E. Garrick and S. I. Rubinow 


SUMMARY 


A connected account is given of tiie Possio theory 
of none tat ionary flow for 3mall disturbances in a t wc- 
dimensiorsal supersonic flow ar.d.of Its. application to 
the determination of the aerodynamic forces on an oscil- 
lating airfoil. Further application is made to the pro- 
blem of wing flutter in the degrees of freedom - torsion, 
bending, and aileron torsion. Numerical tables for 
flutter calculations are provided for various values of 
the Mach number greater than unity. Results for bonding- 
torsion wing flutter are shown in figures and discussed. 
The static instabilities of divergence and aileron 
reversal are examined as is a one-dogroo-af-f roedorn case 
of torsional oscillatory instability, ■ - ■■ 


INTRODUCTION 


The problem of flutter or aerodynamic instability 
for high-speed aircraft is of considerable importance 
and lienee interest is directed to tho aerodynamic problem 
of the oscillating airfoil moving forward at high speed. 
Although for conventional aircraft the subsonic and the 
near-sonic or transonic speed ranges are still of main 
interest, the purely supersonic speed range is becoming 
increasingly significant. 

A theoretical treatment of -the oscillating airfoil, 
of infinite aspect ratio, moving at supersonic speed 
has been given by Possio (reference 1). This treatme'ht 
is based on the theory of small pertmrbations to the 
main stream, thus is essentially an acoustic theory, and 
leads to linearization of the equation satisfied by the 
velocity potential. The airfoil i3 therefore assumed to 
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be very thin, at small angle of attack, and the flow is 
assumed nonviscous, unseparated, and free from strong 
shocks. 

The small-disturbance linearised theory, being much 
less complicated than a more rigorous nonlinear theory, 
is to be regarded as an expedient which allows an initial 
theoretical solution. The theory permits the occurrence 
of weak (infinitesimally small) shocks and .thus the basic 
trends and effects of the parameters of the simplified 
problem can be indicated. The theory reduces to that of 
Ackeret in the stationary (static) case and, like it, 
is not expected to be valid too near M = 1 . In view of 
the restrictions and assumptions in the analysis important 
modifications may be required m certain cases for thick 
finite airfoils, but even here the simple theory for thin 
wing sect! 021s may serve as a basis. 

In addition to Possio's brief work an equi valent 
extended treatment has been given by Borbely (refer- 
ence 2) which utilizes contour Integrations to carry out 
the solution of the partial differential equation for 
the velocity potential according to the Heaviside operator 
method or Laplace transform method. Recently, another 
equivalent treatment has been given in England by Temple 
and Jahn employing the method of character! sties. In 
reference 1 a few curves are given for the aerodynamic 
coefficients but no numerical values are tabulated. 
Reference 2 contains no numerical results. Temple and 
Jahn recognize the lack of numerical results and supply 
some Initial calculations for the functions necessary 
for flutter calculations. 

A paper has recently appeared’ by Schwarz (refer- 
ence 3 ) devoted to computing and tabulating the key 
mathematical functions that arise in the theory. The 
present paper makes use of reference 3 to supply more 
extensive numerical tables for application of tho theory. 
The formulas of the theory are recast: In more familiar 
form for application to the flutter problem and a series 
of calculations on bonding- torsi on flutter are carried 
out and discussed. The performance of 3 imilsr calcula- 
tions for wing-aileron flutter is indicated. Brief 
discussions also are given of the static instabilities, 
divergence and aileron reversal, and of a one- degree- of- 
freedom torsional oscillatory instability. 



NACA TN Ho. 11 5 S 


3 


For complete ness, a connected account of the Possio 
theory ia presented since the original presentation in 
Italian is quite terse and also since it 3.3 believed 
that this treatment Is the simplest ' and most suitable 
for general extensions. The extension of Its applica- 
tion to include the aileron Is given. 


A IT? FORCES AND MOMENTS ON AIT OSCILLATING AIRFOIL MOVING 
AT SUPERSONIC SPEED IN TWO-DIMENSIONAL FLOW 
• Differential Equation for the Velocity Potential 


The differential equation satisfied by^he velocity 
potential In fixed coordinates . in the case of infinitesimal 
disturbances is the wave equation 


JL = ^ 2 0 

c^ b r/ 


(1) 


whe re c 
medium, 
dp 

C*~ = : 

o’p 


is the velocity of sound 3m the undisturbed 
(For the adiabatic equation of state 


Y 


Deferred to a system of rectangular coordinates 
moving forward at a constant supersonic speed v in 
the negative x-directicn the wav? equation satisfied 
by the velocity potential in two-dimensional flow becomes 


JL_ u + £v b Z 0 
c 2 bt 2 c2 cjx bt 




bj 2 


0 


( 2 ) 


It is proposed to treat the effect of aslightly 
cambered thin airfoil moving forv^ard at a supersonic 
sneed v at small (zero) angle of attack a3 that of a 
distribution of small disturbances placed along, the 
x-axis and hence t.o utilize equation (R). ‘Hie velocity 
components in the x- and y-airections relative to 
the moving airfoil are, respectively. 



k 


FAC A T'i l;o. 1153 


and 


b_0 

bx 


'7 


*1 

b 7 


which may be considered the additional components to 
the main stream due to the. disturbance created by the 
oresence of the airfoil. Relative to coordinates fixed 
in space the velocity components are v + v x and v y . 


Effect of a Source ' . - 

Equation (2) 1 b linear and solutions are therefore 
additive, an important particular solution of equa- 
tion (2) having the property of a 'source pulse is 


0 = *-— 
0 V 


A(£, p. T ) 


c 2 (t - T) 2 - - jx - - v(t - T)J 


(3) 


(y - u)‘ 


This solution may he considered to give the effect at a 
point (x, y) at time t of a disturbance of magni- 
tude A originating at a point p) at an earlier 

time T, The potential 0~ o is thus a retarded potential 

and the elapsed time at ( x , y) since the creation of 
the disturbance is T = t - T, ; 

Unlike the situation for a subsonic flow^ Tor a 
supersonic flow the effect of the bis turbance is propa- 
gated only downstream, that is, the point oeing 
influenced (x, y) is always considered to be aft of the 
point of disturbance ( g, Vj). Equation (5) is thus 
valid in the angular region with vertex at ( p) and 
bounded by tv/o straight lino 3 . making the Mach angles 

iu = lsin“^p: = ±sin~ ^ with respebt to the x-axis. 

V .;vJ 

(Gee fig. 1.) Upstream from this angular region tl'jye 
value of 0 O is aero. It follows! also that disturbances 
in the wake need not be considered' and the solution to 
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the boundary problem may be attempted by a distribution 
of potentials of the type 0 taken along the projection 
of the airfoil on the x-axis. 

A disturbance at (£, rj) created at time T is 
first felt at a point (x, y) after a certain time 

has elapsed. The point (x, y) penetrates the wave 
front of the disturbed region and because it is moving 
at a speed greater, than that of the wave front it 
emerges from the disturbed region at a later time . 

Thus, the duration of this initial disturbance at 
(x, y) is r ^ - T . (See fig. ?. *) The transition 

at (x, y) from a region of quiescence to a region of 
disturbance and vice versa is associated with the 
vanishing of the denominator in equation (3)* The 
values of t-j and T 2 for a disturbance created on 

the axis p = 0 are thus given by 

c,) * \j (x - £)2 - y2( K 2 - 1 ) 

; : (4) 

c(M 2 - i; 



where the minus sign Is associated with and the plus 

sign with and where M = <£• It may also be observed 

that a negative quantity tinder the radical sign In 
equation ( 3 ) Is to be interpreted as associated with an 
undisturbed region, (that is, with 0 = 0 ), 


Potential for a Distribution of Sources 

The total effect at any point (x, 7 ) is the sum 
of the effects of disturbances originating between the 
leading edge £ = 0 and the intersection of the Mach 
line through (x, y) with tho Praxis 

S = c-j_ = x - y\/i<i2 - 1 

(since only disturbances created forward of the JTach 
angle region can affect (x, y); see fig. 3 ). 
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The total potential at (x, y) at any time t is 
thus given by ; 


0(x, y, t) = 


i£»l 


nr. 


Jo Ut 


A(£, 0, t - r) 


_ c 2 


Y c2 t 

2 - (X -:s 

& 1 

P T2 A(g, 

1 0 t. 

“1 

dT'i 

1 

' ll 
1 


dr dg, 


dr d£ ( 5 ) 


Boundary Condition and Strength of Distribution 

The function A(£, 0 , t - t ) giving the magnitude 
of the source distribution i3 now to 'be determined by 
the usual boundary condition of tangential flew along 
tho airfoil. If the ordinate of any point” of the moan 
line defining the airfoil is given as y = y m (x, t) the 

boundary condition may be written 


— w(x, t) = 

\^Vy=0 d ' t 

bvrn : 

ox 6 1 


whore w(x, t) thus represents the vertical velocity 
induced by the source distribution in order to realize 
tangential flow at the airfoil boundary. (In the 
stationary case - Ache ret treatment - the two surfaces 
of the airfoil may be considered ais acting independently, 
which can also be done for the nonstationary case. 
However, for the purpose of obtaining the oscillating 
forces in the linear treatment it Is sufficient to - con- 
sider separately tho upper and lower sides of only the 
mean line. ) 
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The evaluation of -rf- as y approaches zero may 

°y 

be readily obtained by use of the variable ' 0 Instead 
of t where 2 t ~ (t£ - T]_ ) cos 0 + + T-j_. This 

substitution in equation ( 5 ) yields 




T 2“ T 1 


cos 0)d0 d£, 


By differentiation with regard to y and with the aid of 
an integration by parts - -"j 


9l = ^ 

uy y'v 2 - - c 2 5y 



0, t 


My \ 

CV& 2 - 1/ 


+ 


y 


■y/jrZ _ C 2 C -y/f V i2 _ 1 



^4 sin 2 e d0 dc, 
6t 2 


Since ~ x “ y\-4l^ - 1, there results in the limit 

as y approaches zero on the positive side, the important 
relation 


W\ 


IT 


, ~ J ~ A(x, 0, t) 


y=+o 


or, briefly. 


A(x, t) =- — w(x, t) 


( 7 ) 


Por y approaching 0 on the negative side .an 
equal and opposite result is Detained and hence the 
distribution of singularities to be utilized to replace 
the airfoil is of the source-sink type. Thus $ is to 
be understood in the subsequent analysis to be prefixed 
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by a - sign, + for the upper aide and - for .the 
lower side* 1 


The total potential for y = 0 may now be expressed 
by means of equations (5) and (7) as 


0 (x, 


1 1 j )x P T 2 w(g, t - T ) 

W - 1 Jo Jt 1 - Ti)(t 2 - T j 


dT dg (8) 


where, from equation (I4.) vith y — 0, 

v - g 1 : 

T 1 - — * *- 

C M + :1 


and 




Application to Oscillating Airfoil 

The general result given by .equation (8) may now 
be applied for. definiteness to the: case of an airfoil 
performing small sinusoidal oscillations in several 
degrees of freedom., let the wing iundergo the following 
motions: a motion due to displacement h (velocity h) 

in a vertical direction; a torsional motion consisting 
of a turning about x = xq with instantaneous angle of 

attack a; a rotation of an aileron about its hinge at - 
x = x-j_ with instantaneous aileron angle p measured 

with respect to a. (See fig. k.), 

In accordance with equation (6) the vertical velocity 
at any point x of the airfoil situated at 0 ^ x ^ 2b 

(of chord 2b and leading edge at |x = 0) is easily 
recognized to be ; 

w(x, t) = - |k + va + (x - x 0 ^a + vp + (x - (9) 

where the p-terms are to be Interpreted as zero for 
x < x^ (and where the minus sign ;is introduced because 
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the vertical velocity w is positive upwards whereas 
the terras within the brackets are positive downwards). 

It is convenient In treating sinusoidal motion to 
utilize the Comdex notation 


h 

a 

P 


= iiQe 
- CCQS 

= Po e 


ICOt 


icot 

icot 


( 10 ) 


where hQ, clq, and {3g are complex amplitudes and 
hence include phase angles. 


Since the further analysis is concerned only with 
exponential time variations of the type given In equa- 
tion (10), the function w{ g, t - t) occurring in 

equation (8) is of the form w(g).e ia> ( t "' T ^ , which' may 

also be written for convenience as w ( g, t)e -iwT . The 
potential 0 given by equation (8) may now be written as 


j#(x, t ) = - 


\M Z - 


\i 




w(£, t) I(£, x) dg (11) 


0 


where 


!(£, x) 



The integration with regard to t may be readily per- 
formed by substitution of the variable 6 where 
2 t = (t£ - Tq) cos 0 + T£ + t Then. 


m 

I (^ x ) = I e -iw(T £ +T 1 )/2 / ^-icocosGC^-T^J/a 


d8 


to 
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With, t-j_ and T 2 replaced by their values as given for 

equation ( 8 ) and with the aid of the Bossel function 
relation > 


1 

TT 



e -iXcos 0 


d 0 


= J 0 (\) 


It is recognized bhat 


, x-g' m / „ 

-iCO — - — — (X — ci 

Kg, x) = e ° M2-1 J Q [- S 


0) 

M 2 - 1, 


( 12 ) 


Throughout the subsequent analysis it is convenient 
to employ the variables x and jg In a new sense to 
mean nondimensional quantities obtained by dividing the 
old variables by the chord 2 b. Hie retaining of the 
symbols x and g for the nondimensional variables 
shoixld lead to no confusion. ! 

The potential 0 of equation "(11) is then 


0(x, t ) = 



jva + h + 2 b (g - Xq^cx + vp + 2 b (c, 




where with the introduction of the important frequency 
parameters ; 

, _ ; 
v 


(0 - 


2h-L 


M £ 


~1 


the function I(g, x) becomes 

I(g, x) = e --to(x-g) j 


0 


:co 

M 


(x - g) 


( 12 1 ) 


C 
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Thus, I( £,, x) is a function of the variable x - £, and 
of two parameters M and to, or, alternatively, M 
and k , 

It is desirable to express the potential 0 as the 
sum of the separate effects due to position and motion of 
the airfoil associated with the individual terms in 
equation (13). Thus 


0{x, t) = 0 a + + 0'a + 0$ + 


04 ) 


where 


2b 


*tx 


0a = —, nz: - ~r . va I(g, x) d£ 

V& 2 - 1 Jo 


_ 2 b_ 

Vwl 2 - 1 




^h = 11 I dg 


/O 


4 = 




\4 2 - 1 


a 


(g - * 0 ) Kg, x) d£ 


o-h p x 

yfi i(£, x) d£ 


K - 1 Jx. 


0k = “= pf (S * * 1 } KS, x) d£ 

vk 2 - 1 ,JL 
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Forces and Moments 


The basic nressure formula in; the theory of small 
disturbances is 


P = 



which in the present case of the moving airfoil may be 
expressed as ; 


P = 



where p is the density in the unjdisturbed medium. The 
local pressure difference on the airfoil surface betwoon 
the upper and lower surfaces at any point x ( nondimens ional) 
is 


P' 


-2f 


60 v 

6 1 2b dk> 


(15) 


The total .force (positive downward) on the airfoil is 


P 1 

p - 2b 1 p 1 dx 

Jo 



r-i; ' 

Ippb I 0 dx 

Jo I 


(16) 


The moment (positive clockwise; fig. ip) on the 
airfoil about any point x Q is ■ 



entire 


( 17 ) 
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Similarly, the moment (positive clockwise; fig. )(.) 
on the aileron about the hinge point x^ is 



In the further reduction of equations (l6) to (l8), 
with the potential 0 replaced by its separated form 
given in equation (1 /l), the following seta of integral 
evaluations are required: 



60 , 


a . 


dx 


dx 


2b 


y£i 2 - 1 


vor-j_(M, 


k) 





k. 





2b 


J a 


dx = — =r 


yfl 2 _ 1 


vargCM, k) 


0a dx 


i4b^_ 
\/il 2 - 1 



k) - x Q r 2 (M, k). 
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a 

dx-i 


dx = vp t2 ^- k - * 1 ) 


2f* dx. 


- _J±^ 


- P § t : $( M » i k > x l) 


v^T K 2 * 


*> 0 a „ 2b ,„ , 

x dx = — — — — ^ va q (M, k) 

6x y^ 2 - 1 - 1 




a 


x dx 


■ jjhj 


/o dx y^ 2 - 1 


r 1 

c. x dx 


r= a 


| q 2 (M, W) - x Q q 1 (M, k)J 


Ux x ox 


a . 

— c - x dx 




_ lib 4 


====== V (\ j~3 1 (^ k, Xy) +-X X t x (»!, k, X 1 )J 


/x l 


dx yM 2 - 1 


I s 2 ( M » k > *h) + *1 t 2 (m, k, x 3 ) 


r 


2b 


0 a x dx = —~== va w q ? (M, k) 
yM^- - I - ; 


0.x dx 

r nr 


a 


I|.b 2 . 

a 

\L 2 - 1 


1 , x a 

~ q ? (M, k) ~ 2 x O ( 1 2 (M ’ k) 


I 
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dx 








P 5 ( r > k >*l) - |x 0 p z( K ' k ’ x l)j 

\ k ' x l) 



The functions defined by the {foregoing Integral 
evaluations are further discussed i in the following sec- 
tion; first, however, the force arid moments (equations (16) 
to (18)) are given in their final ! forms as 


- _ 


yjy[2 - 1 L 


v^va + h - 2bxQ&^r-j_ + - 2b (iva + h; - 2bxQa^r2 


o r 3 T 

+ Iqb^a — - + v'-pt, + l+bvpt? + ]. |_b^p ~ 


( 16 * ) 


M a =■ 


80 b 2 


v(va + h - 2bx 0 a^)qq + 2b(dva -t h - 2bxQ<i^ 


fivl 2 - 1 L 

q, ^ p 

+ Ub 2 a yf- + v 2 p(s-, r x.-, tq)’ ■«• l.|bvpy~ t x^ih, 

/ * “ 

+ l+b 2 p(^|. + xi y 2 ) - 2bx d P 


( 17 ' ) 
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M. = - ~~=rj~v^va + h - 2bXQcJp^ + 2b^2va + h - 2bXQaJ — ^ - 
^ ym^ - 1 L 0 


+ Iqb^a + v^p^ + 1+bvp — + l).b p -g 


2-* 3 3 

% a ..i * ■ 


2 

(18*) 


Reduction and Evaluation of Poregoing Integrals 

It is convenient to introduce the substitution 
u = x - £, and to express the function 1(5,, x) 
(equation (12 1 )) as 


/— 


I(i, x) = I ( u ) = e“ la)u J 0 (| 


u 


(19) 


The various functions defined by the foregoing sets of 
integrals may now be expressed as follows; 

I' 1 

r^M, k) = I(u) du 

Jo 

nl 

r^(M, k) = I I I(u) du dx 

Jo Jo 

f r- 

r^(M, k) = 21 I (x - u) I(u) du dx 

J) Jo 


nl 


q^(M, k) = j u I(u) du 
u'O 

nl j* 

q_ ( M , k ) = 2 j I x I ( u ) du dx 

Jo Jo 


P 1 

JO (JO 


q^(M, k) = 6 j j x(x - u) I(u) du dx 
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k > x i) = | (u - I(u) du 


ux n 


P 2 ( m , k, x-J = 


f')i 0 * 

2 (x - x^) I(u) du dx 



f)l nx 

P 5 (m, k, x.J = 6 (x - x x )(x - u) l(u) du dx 


UXj _ (JO 


i-xi 


k, x-jJ = I I(u) du 


JO 


^2 x-^y — 


^(M, k, xj = 


ii-x-l nx 


I(u) du dx 


0 JO 
'>1-X]_ px 


(x - ;u) I(u) du dx 


Jo 


1-xi 


k, x-jJ = u I(u) du 


s 2 ( m , k, x x ) = 


Jo 

pl-X]_ px 

2 / / x I(u) du dx 


JO JO 


3 (m, k, x-J = 6 


x(x - u) L(u) du dx 


JO 
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Borbely (reference ’2) has shown by' means of 
reduction formulas that the six r- and q-functions may 
be obtained from a single integral. . In a similar man ner 
it may be indicated how the foregoing 15 functions may 
be obtained from the evaluation of the same integral. 

The reduction Is accomplished in two stages. First, 
consider integrals of the following type : 


- fy(M, co) = I I(u) du 
>0 


S\ = 


(m, (oxj = 


X+l 


I(u) u^ du 


> ( 20 ) 


“1 


h X “ f \ 


M, u(l - x-J] = 


i 1 " *l) 


X+l 


Al-X- 


10 


I ( u ) u^ du 


By Integration by parts It can be readily verified that 
the following relations hold 


r l = f 0 

r 2 = f 0 ~ f l 

r 3 = f 0 " 2f l + f 2 



q^ — 2 Tq - 3fq + fj 
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p i = q i " x i r i + x i 2 (so “ Si) 

p 2 = q 2 - 2x ± r 2 + X 1 5(g 0 - 2g 1 + g 2 ) 


q 3 • 

- 3 x irj 

* x iK«o 

- 3 Sx + 3 g 2 

v 

= (- 

x l ) h 0 

• i 

t 2 

= (l - 

x l ) 2 ( h 0 - 

h l) 

fe 3 

= ( X ‘ 

x l ) 5 ( h 0 - 

2 h 2 + h 2) 

a l 

= (1 - 

x l) 2h l 

■ s ; 

s 2 

= (i - 

x i) 3 ( ll o - 

h 2 ) 

s 

3 

= i 1 ~ 

x l ) l! ( 2h 0 

" + hj) 


The final stage In the reduction of these functions 
is to utilize the following recursion formula (refer- 
ence 2 ) obtained by integration by parts: 


M - ] _ , 

— to) = 

M" 


1 + (1 - X) -~ 

to 


e" 1W j r± 
0\M, 


1 « - 1 CO (M 

5 e j iU 


+ i ( 1 - 2\) f T, ■ , ( X, CO) 
' '< ~ .1 . 


+ (1 - X) 2 ~ f’ ( III, co^ 
to X-2 


( 21 ) 


where X. ^.1 and f with a negative subscript is to be 
interpreted as zero. - ! 
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The function f\(M, 03 ) may clearly refer also to 
the foregoing g- and la-functions, if cB Is replaced 
by the appropriate parameter; namely, cjox q for 

and «(1 - s 0 ) for . (See equations (20).) The 
recursion relation (equation (21)) thus reduces the 
various functions to the single function 





co 


UO 


e -iu 



du 


( 22 ) 


which is therefore the only integral needed in the 
evaluation of the forces and moments. 


The important Integral in equation (22) has been 
recently made the subject of a mathematical investiga- 
tion by Schwarz (reference 4) . Schwarz gives tables 
of the values of its reed and Imaginary parts to eight 
decimal places for 0 = oc £ 5 and for 1 ^ ^ 10 for 

conveniently small intervals. For values of To > 5 not 
given in Schwarz' tables, the function fg may be 
evaluated by means of the following series development 
(reference 2): 




1 

l(2n+.l) 



+ iJ 


n+1 



selected values of the Mach number 

o 5. 

J z’__ 

of 0 ) 


M = — - 


Table I give 3 values of the functions f c (M, co) 

based on the tables of Schwarz and on equation (23) for 

10 

10 9’V. 7’3' 

- — y and 5 and for various appropriate values 

3 / l\ 

ior y ). Later use is made of the values given 
in table I for obtaining tables for flutter calculations 
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EQUATIONS OP MOTION AND DETERMINANTAL EQUATION FOR 

FLUTTER CONDITION 


The equations of motion and the border-line condi- 
tion of unstable equilibrium yielding the flutter speed 
and frequency may be obtained exactly as in the incom- 
pressible case treated, for example, in reference I 4 ., 

The two-dimensional treatment (infinite aspect ratio) 
is retained herein. Modifications: due to assumed vibra- 
tion modes of the finite wing may ;of course be introduced 
a3 in current practice (for example, reference 5). The 
modification of the forces and mcmtents due to the three- 
dimensional nature of the flo\7 is a more difficult 
problem which remains to be studied. 

The equilibrium of the vertical forces, of the 
moments about the torsional axils x = Xg, and of the 

moments on the aileron about its hinge x = yields 

the three equations, > 




hM + 

ds a + ps p + hC h = P ! . 

. i 

! 



*dl a + 

’pjlp +- 2 b(x 1 - x 0 )3p] + h3 a H- 

aC a = M a 

>- 

(2k) 

Pip + 

a Jlp + 2 b^X]_ - Xg^Gpj + hSp .+ 

pc p = Mp 




where the various parameters arc defined in the list of 
notation. (See appendix.) 

In order to define the border!- line condition of 
unstable equilibrium separating damped and undamped 
oscillations, the variables h, or, and p are used 
in the sinusoidal exponential form, given In equation (10), 
For the desired condition, it is hocessary that the 
equations ( 2 ) 4 .) have a (nontrivial) solution for the 
complex amplitudes hg, fig, and ; Pg, or that the 

foT lowing detorminantal equation hold: 
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A ch 

A ca 

A c p 



A ah 

Aaa 

A ap 

= 0 

( 25 ) 

A bh 

^ba 

Abp 




where the complex elements of the determinant In 
separated form are 

A ch = Q&x - w + L X + 1L 2 

A ca = -V*a + l 3 + iL [}. 

A cp = -tap + 1*5 + iL 6 

A ah = ‘^ x a + M 1 + m 2 

\a = Q aX " M-r a 2 + H3 + 

A a|3 = (fp 2 + 2(x l “ x 0) x p] + M 5 + “6 
Ab h - -IJXp + % + 1W 2 

Aba = [ r p 2 + 2(x l. " x O )x p] + n 3 + 1N U- 

A bp — + ^5 + -^6 ’ 

and where the L's, M's, and IPs are defined by the 
force and moment equations (l6 ! ), ( 17 r )? and (l8 l ) 
expressed In the following forms: 
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P = -i^pb v 2 k 2 e lWt }(^ (Li + ilg) + a 0 (l^ +-1%) + P 0 ^+iL 6 )| 

M a = -kpb 2 v 2 k 2 e iWt (~)(Mi + iM 2 ) + + P 0 (j% + iM^)j >(2 6) 


Mp =. -ij.pb 2 v 2 k 2 e i6)t 


Hence , 



*1 + !%) * «o(b i + iN i,) + ft>(% + l *c) 


L 1 + 1L 2 - ,3—- {-2-2 + | 


Lj + ILj^ — 


V^ 2 - 1 
1 

.'Ivt 2 - 1 


o . 21 

-2r T + — 

3 li 


r i) . 

r r a - s( 2r 2 + f r i) 




-2x 0 \-2r 2 + - r. 


^ + 1L 6 “ 


•il 2 - 1 


-at, + r= t, - |(-2t 2 + I t x 


21 

Jj T '* U-T 

3k 2 


M x + iM 2 
m 3 + iM l|. 


\il 2 - 1 
\4 2 - 1 , 


21 

l 2 + k 
21 


~ 2q 2 + T q l ” 2x 0\“ 2r 2 + 

1 

1 / 21 


k r l 


"5 q 3 + T q 2 “ S lr 2<1 2 + T “i. 

21 l/_: i 


21 


_2x 0 I“ 2r 5 + k r 2 "-,k \ Zr 2 + k 1 !)” 2 ^ + k q l 


~ 2 x 0 (- 2 r 2 + 


* 0t 


M r + iMk 

p D 


Ml, 


21 


+ * 7 — 3p — rrl~2'3 0 + 

/k2-i 3 3 k 2 k V 2 k 1/ 


21. 


: ( x i- x o) j- 


+2(Xi -Xo) >- 2t^ + — t 

k 




L.\il 


2 - |Er 2fc 2 + i 
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% + iN 2 
N 5 + 1Nj,_ 




N 5 + iN 6 




The de terra inantal equation {Z'y) with the foregoing 
complex elements is equivalent to two real simultaneous 
equations and hence may be solved for two unknowns . In 
a given case the usual unknowns are' the flutter speed _v. 
and the flutter frequency w or, more conveniently, the 
related ncndimensi anal parameters X and l/k. The 
paramoter X apoears linearly and only in the major 
diagonal elements (with bars), while the parameter l/k 
appears transcenden tally in every element of the deter- 
minant. Hence an obvious procedure though not the 
simplest for obtaining the simultaneous solutions of 
the two equations is to fix values of l/k, to solve 
for the roots of the two polynomials in X, to plot 
graphically these roots against l/k, and to note the 
points of intersection. 

In a systematic numerical study of flutter any two 
parameters maybe utilized as unknowns instead of X 
and l/k, which is often more convenient. A discussion 
of such procedure and the use of a method of elimination, 
for simplifying the calculations is given in the 
appendix of reference 6. 

The application to tho two-degree- of -freedom subcase 
of bending-torsion flutter is treated more fully in the 
following section. 
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APPLICATION TO BENDING -TORS ION FLUTTER 


The determinantal equation in the two degrees of 
freedom h and a is 


■oh 

A oa 

ah 

A aa 


0 


or 

— (x + L z + iL 2 
-l^a + M 1 + 1M 2 


“l^- a + + iL^ 

0 a x - txr a a + + iMJj 


= 0 (27) 


The two equations in X obtained by equating the real 
and imaginary parts separately to:zero are 


%^ a X 2 + j^ a (Li - q) 

and 



(pa l 2 + + Cj = 0 


0 


? (27') 


where 

°R = t J 'fa( M l' + l 3) ~ ( m 3 " ‘ Ir a 2 ) - %^ a 2 



+-d r 


and 



L 2^ 2 


+ 



where 


D r - LqMj - L^M^ - L^Mj^ + L|^M 2 


D I = - L 1( _M Z + L £ M 5 - L 5 M 2 


and 
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For convenience in numerical tabulation, it is 
desirable to introduce primed quantities, independent 
of the parameter Xg, defined by the following relations: 


L 3 

N"\ 

II 

2x Q I 

1 



L[(. 

II 

- 2 XQ 1 

'2 



Mi 

= My 

- ^o 1 

'1 



m 2 

= My 

- 2x 0 I 

'2 



M 3 

= V 

l 

£ 

0 

[ ( V 

+ iy ) 

“ 2x 0 L l 


sr 1- 

II 

" 2*0 


+ 

" 2 * 0 L 2 


(28) 


In table II convenient expressions for the quantities 
L x , L 2 , iy, y, M 1 », My, My, and My are 

given and tabulated together with the c cmbinations 
M]_* + iy and My + iy . Clearly these quantities 

depend on the function f Q given in table I and hence 

the tabulation is made for the same values of M 
and l/k (or 5J) , In addition, table II contains 
values for the quantities Dr and Dj which, in fact, 

are independent of x Q and may be expressed as 

d r = ryy - ry M y - Lyy + 

and 

d-j- = iy/y - + iyiy " l 3 ,m 2 t 


The numerical application in the case of bending- 
torsion flutter has been performed for .various selected 
examples. In most of the calculations the numerical 
procedure was to fix values of l/k, eliminate X, and 
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solve for the parameter x a . Interpolation was also 

used to obtain additional points In order to Improve the 
fairing of some of the curves. Values of l/k less 
than 1 did not yield any flutter points in this 
procedure. Results are shown plotted in a number of 
figures (figs, 5 to 20 ); however, before these figures 
are discussed, it Is .desirable to explain the significance 
of the parameters and the numerical values assigned to 
them , 


The parameter p may be considered to signify the 
wing density and three selected ; values 3*927* 7 • 85*4-* 
and 15,708 in the order of increasing wing density have 
been mainly used, in the calculations, (These values 

correspond to values of .-j = 5 * 10 , and 20 in the 

notation of reference J4. , } Alternatively, an increase 
in p may be interpreted as an 'increase in altitude for 
a fixed wing density. The parameter' p may be expected 
to range up to high values for actual supersonic wings 
at high altitude. Only a few calculations, however, 
have been made for high values of p ( p = 78*5^-* 

- “ 100; see fig. '18) , ! * 

The parameter 00 /00 is thje ratio of the wing 

bending frequency to the wing torsional frequency and 
may be expected normally to be less than unity. The 
throe values 0, 0,707, and 1 havte been largely used in 
the calculations although other .values up to 2 have 
also been studied. 

The parameter x Q represents the position of the 

elastic axis measured from the leading edge and the three 
values 0,1}., 0,5, and 0,6 represent, respectively, posi- 
tions at lj. 0 , 50 , and 60 percent chord, (These values 
correspond to values of a = “0,2, 0, and 0.2 in the 
notation of reference I4. , ) 

The parameter x a represents the distance of the 

center of gravity from the elastic axis. For example, 
x~ - 0,2 represents a position of the center of gravity 
10 percent of the chord behind the elastic axis. In 
many of the calculations x a has been regarded- as 
variable. ■ - ; 
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The parameter represents the radius of gyration 

of the wing about the elastic axis and has been kept 
fixed at the value r a 2 = 0,25. 

The ordinate in figures 5 to 20 is the nondimens ional 
flutter coefficient v/bm a where ba) a is a convenient 

reference speed. This coefficient is also a function 

of the Mach number M = — .and several values of M 

c 

have been employed In the calculations 

In a plot of the flutter coefficient v^bo^ 

against M, straight lines drawn from the origin at 

angle 6 and Intersecting the curves may be given an 

.interesting interpretation (fig, 17 ) * The slope of " the 

v/bo j o bco a 

line Is given by = — - or cot 5 = . Thus, 

v/c bco a c • 

cat 6 is directly proportional to the product of the 
chord and the' torsional frequency. The question of 
whether at a given value of M the value of bo> a 

which wi 3 „l just prevent flutter is also sufficient to 
prevent flutter at neighboring higher values of M is 
answered by the simple criterion of whether cot 6 
Increases or decreases. In figure 17 two typical 
flutter curves are shown. In curve B the value of bco a 

just necessary to prevent flutter at a speed corresponding 
to the value of M at P£ is insufficient to prevent 

flutter at any higher value of M for which the flutter 
curve is below the straight line OP^, For the type of 

curve A a maximum value of S occurs at the "design 
critical points" P^. The values of boo a just necessary 

to prevent flutter at a speed corresponding to the 

value of M at P-> is also sufficient to orevent flutter 

at all higher speeds. 

The following salient facts may be extracted by 
Inspection of the figures. Flutter exists or is possible 
for various ranges of the parameters but, in general, 
compared with subsonic cases the ranges of the parameters 
yielding flutter are more restricted. 

The chordwisa position of the aerodynamic center, 
the center of the oscillating pressure, is an important 
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factor in the consideration of flutter* In the simple 
theory the mid chord is the aerodynamic center for 
M » 1. For subsonic speeds, M « 1, the linearized 
theory indicates the quarter-chord position as the aero- 
dynamic center. It should be expected that in the 
transonic region near M = 1 the aerodynamic center 
may shift considerably. From this point of view alone 
conclusions drawn from the simple theory for the range 
near M = 1 may require large modifications. The 
nature of the modifications may be roughly Inferred by 
further experimental and theoretical study of the 
behavior of center-of -pressure locations. 


For low values of the ratio of bending frequency to 

Wh 

torsional frequency 77 - ~ 0 the position of the center 

of gravity relative to the aerodynamic center is Important, 
For center-of -gravity positions forward of the midohord 
no flutter exists, whereas for positions behind the raid- 
chord there Is a sharp decrease in the flutter coeffi- 
cient from infinity; the position of the elastic axis 
influences the value of the flutter coefficient in this 
region, forward positions being more favorable (figs, 5 (a) 
to 16 (a)). 


For values of 



the position of the center 


of gravity relative to the elastic axis becomes of more 
Importance. For center-of-gravity positions forward of 
the elastic axis no flutter exists, whereas for positions 
behind the elastic axis flutter does ocour, and a relative 
minimum coefficient appears for center-of-gravity posi- 
tions only slightly (a few percent of the chord) behind 
the elastic axis. 


The intermediate case, for which . — - = 0*707, 

w a 

shows a blending of the effects in which the center- 
of-gravity position relative both. to the aerodynamic 
center and to the elastic axis is significant. 


In figures 12 and li^. there are shown, for refer- 
ence, some numerical values of- co/co^, the ratio of the 

flutter frequency to the torsional frequency. 
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The effect of the wing-density parameter jjl is 
rather complicated but, in general, an increase in p 
yields a corresponding increase in the flutter coeffi- 
cient. For low values of co^/co^ and for high wing 

densities this increase is expected to be proportional 

/— w h 

to I 11 the resonance -like region near — = 1 

& 

and for small values of x a the flutter coefficient is 

relatively unaffected by the value of n, and in this 
region the structural damping may be expected to be 
particularly effective in increasing the flutter 
coefficient . 


For values of the Mach number near unity (for 
example M = a value for which the validity of the 

theory is in question), the flutter calculations become 
difficult to plot because of the appearance of other 
branches. In some cases (for instance, x n = 0.6) the 


flutter instability appears limited to a definite range 
of speed. Calculations to include damping were performed 
to verify the existance of the range. “(The appearance 
of these other branches seems to involve values of l/k 
for which the quantity Mj^ is negative. The condition 


of negative is significant for the one -degree- of - 

freedom instability discussed in the next section. 


A plot of the flutter coefficient against Mach 
number for two values of x a is shown in figure 17. 

The significance of the straight lines drawn from the 
origin has already been discussed. The type of curve' A 
is representative of the effect of forward location of 
the center of gravity and the type of curve B is 
representative of rearward locations of the center of 
gravity. Figure 18 gives a plot of the flutter coeffi- 
cient against M for various values of the wing-dehsihy 
parameter \i and for a rearward location of the center 
of gravity. The subsonic values for M = 0 and 
M = 0.7 shown on these curves _and on some of the other 
figures have been either taken from reference 7 or 
calculated in the manner outlined therein. The subsonic 
and supersonic parts of the curves (figs. 17 and l8) 
have been arbitrarily joined by a dashed smooth curve In 
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the transonic range. In figure 19 there is given a 
cross olot of flutter coefficient against frequency 
ratio ( 0 ^^, for various values of M, and in fig- 
ure 20 is given a similar croas plot for three values 
of the elastic- axis parameter x Q . 

An indication of the effect of structural damping 
in Increasing the flutter speed in a few examples nay 
bo obtained from the following table, where g a 

and g. are the torsional and flexural damning coef- 

10 

ficienta, respectively, and where M = — , ^ = 7 

a = 0, and x a = 0.2: 


co h/ co a 

, 

8h 

"/“a 

v/bco a 

0 

0 

0 

0.673 - 

2.1+38 

0 

•05 

0 

• 6l+3 

2.551 

0 

.10 

0 

.628 . 

2.669 

»?°7 

0 

0 

• 777 

1-555 

*707 

.05 

0 

• 771 

1.553 

.707 

.10 

0 

.7 b6 

1.569 

.707 

0 

.05 

.788 

l.5?2 

.707 

0 

.10 

•707 
• i / l 

1 . 6l\Z 

.707 

.05 

.05 

.782 

1 . 623 

.707 

.10 

.10 

.781+ 

1.725 


STATIC CASES - '.VI'txG DIVERGENCE AND AILERON REVERSAL 


It Is of some interest to examine the expressions for 
tho forces and moments in the limit case In which the 
frequency anproach.es zero. Thors follow then for the 
mean-line wing section the well-known static-caso results 
which nay of course be obtained directly without the 
use of a limiting process, as originally treated by 
Aekeret„ Thus, with the use of the following relation 
easily verified from equations (20), 

lim f\(m, k) - — 
k. — * 0 \ + 1 
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there are obtained from equations ( 16 1 ) to (l8 r ) for 
the lift and moments in the static case. 



22 


“o =- 


Lj-pb^v 


yfe 2 - 1 *- 


h- 


2xq ) cc + 


(l - xj (l + x ± 



Mr. 


= (l - + 

yki' - 1 


These relations for the mean-line wing section are 
now used to obtain the critioal speeds for the” static" 
instabilities -wing divergence and wing-aileron reversal 
(for wing of infinite span). At the wing’ divergence 
speed the effective torsional stiffness of the wing 
vanishes, hence the total moment about the elastic axis 
is zero. The sum of the structural restoring moment and 
the aerodynamic twisting moment is 


ippb^v^ /, _ \ 

+ a(l - 2x 0 ) 


aC 


a 


i/d 2 - 1 


which when equated to zero yields the divergence speed 

V D = '“aC* 2 - 1 ) 1/ ’ 4 ' VK 2 '• -Ur - 

\.iz 0 - 1 

Thus, the divergence speed is real only for positions 
of the elastic axis behind the aerodynamic center (mid- 
chord, in the simple theory). This formula obviously 
should not be used for values of II too. near unity. 

Per comparison it is of interest to note the 
corresponding result for the divergence speed in the 
subsonic case, where the aerodynamic center is (approxi- 
mately) at the quarter- chord point. Thus, 



ITACA TIT IT o. 1153 ' 




v. 


D 


^(1 - 


\4jCq - 1 


whore M < about 0,7. ■ 

The aileron reversal speed Is deterrmined by the 
condition that the change in angile of attack due tor-wing 
torsion nullifies the effect of movement of tho aileron 
so as to yield no change in lift; (in rolling moment, in 
tho case of finite wing span). There are two equations 
to be satisfied for this condition: namely, 

gc + ^1 - p — 0 
(that is, L = 0) arxd 

aC a + 


i(i - 2x o) a + (1 - x i) 0- 


+ 


- lx 


0, 




Tho aileron reversal speed, obtained by elimination 
of a and p, is 

V H = t® a (i;2 - l) lA vG? ~ 


For _hfnge positions aft of tho midchord, the factor 
l/vxp in this expression varies; from 1 .Ip to 1.0. The 

aileron reversal speod is thus relatively unaffected by 
the nosition of the hinge. In general Vr? may be expected 
to be ]. ovrer than v-n>. : 


0 IMF - DEGREE - OF- FRE'SDO M OSCILLATORY INSTABILITY 


As was pointed out by ?o33iD, the theory indicates 
the existence of a torsional instability which may arise 
for a wing having only one degree of freedom. This 
instability is due to the wing being negatively damped 
in torsion and is associated with the vanishing (and 
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change 'in sign) of the torsional damping coefficient Mj, 

(equation (26)). " 


Certain considerations for the ease -of slow oscilla- 
tiona made, by Posaio (reference 1) and further diseuSSod by 
To mole and John serve to- bring 'out the .main results. 

Thus from equation (20)/ for slow oscillations. 


f x (M, k) 


1 

\ + 1 


. 2kM 2 1 

x 

- 1 \ + 2 


and 


\ = 


1~2 


#T 


1 k 1 3 


k ~ 9*q + 6x 0 2 




The condition XqJ = 0 is shown plotted in fig- 

ure 21, where the shaded area is. the region in which the 
instability is possible (negative M[j_). The maximum 

ranges for the parameters x 0 and M in this region 
are x Q less than 2/3 and M less than y^» 5 (and 
greater than unity).- 


(It may be appropriate to mention that a similar 
torsional instability is theoretically Indicated even 
in the subsonic (incompressible) case for positions of 
the axis of rotation between the leading edge and the 
quarter-chord point. However, the combination of 
parameters required for this indicated instability is 
practically unattainable . ) 


The torsional instability may be studied more 
fully in the general case. It is found that the range 
of instability for the parameters x Q and M remains 

essentially as in the simple case (large l/k) but more 
information may be obtained regarding the critical 
3need and frequency. The moment equation is equivalent 

to A =0, or to the two equations 
aa 

Q a X - pr a 2 + M 5 (m, x q ) = 0 
(m ) x q) + SqPq/- = 0 
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where the structural damping coefficient in torsion g a 
has been introduced as in reference 6. The critical 
speed and frequency may be studied as functions of the 
parameters xq, M* g a and the product combination (j,r a 2 . 

Results of a few selected calculations are shown plotted 
in figure 22, Since instabilities are Indicated for the 
range of near-sonic values (1 < M < I.58) it would 
seem that a more comprehensive investigation of this 
problem is very desirable. 

It may be remarked that a similar analysis for pure 
bending exhibits no instabilj ty while the case of the 
aileron alone does exhibit a range where such Instability 
may occur. This range for an aileron hinged at Its leading 

edge Is 1 < M £ 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , May 29* 
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APPENDIX 

SYMBOLS 

0 disturbance velocity potential 

t time at which disturbance influence is felt 

T time at which disturbance is created 

T = t - T 
p pressure 

p 1 pressure difference 

p density 

Y adiabatic index (for air, y ~ 1 -I 4 . ) 

v velocity of main stream (supersonic) 

c velocity of sound in undisturbed medium 

M Mach number (v/c) 

x coordinate measured in direction of main 3 troam 

y ordinate 

x 0 abscissa of axis of rotation of wing section 

(elastic axis) 

abscissa of aileron hinge 

c, T) abscissa and ordinate of point of disturbance 
b one -half chord 

After equation (12) the quantities x, y, x 0 , xq , 
and £, are employed nondimen sionally and are referred 
to the chord 2 b as reference length. 

w(x, t) vertical velocity at position x on chord and 
at time t 



-r n 
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h vertical displacement of jaxis of rotation 

a angular displacement about axis of rotation 


angular displacement of- aileron; measured with 
respect to a : 


(0 

k 

co 


angular frequency of oscillation 

reduced frequency (cob/v) 

- , / 2k m2 

frequency parameter 


- 


I(£, x) function given in equations ( 12 ) and ( 12 * ) 


J n (l ) Bessel .function, of order In . - 

The following additional symbols , employed in the 
flutter equations, conform to the notation used in refer- 
ences iq and 6 , in which the half-chord b is the unit 
reference length. 


M mass of wing per unit span 

S a static moment of wing-aileron combination per 

unit span referred to the elastic axis 

S3 static moment of aileron per unit span referred 

to aileron hinge : 


I a moment of inertia of wing?- aileron combination 

about elastic axis per; unit span 



moment of inertia of aileron about its hinge per 
unit span 


a 


c 


a 


coordinate .of elastic axis measured from mid- 
chord (2 xq - l) 

coordinate of aileron hinge axis measured from 
the midchord -^2xj - 1) 


location of center of grairlty of wing-aileron 
system measured from elastic axis S ; 

location of center of gravity in percent 
total chord measured from leading edge 


is 100 


1 +- a + x, 


a 


= 100 (*0 


f) 
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a 


'a 


to 


a 


<0 , 


co 


h 


v = r “ 


TT 1 
L K 


reduced location of center of* gravity of aileron 
referred to c (Sp/M-A 

radius of gyration of ( wing -aileron combination 
referred to a 



reduced radiu s o f, gyration of aileron referred 
to c 


per unit 



torsional stiffness of wing around a 
span 


torsional stiffness of aileron system around c 
per unit span 

stiffness of wing in deflection , 

natural angular frequency of torsional vib rations about 
elastic axis (“a = 2 Trf a» where f r 

is in cycles per second^) 


-a 


natural angular frequency of torsional vibrations 

■ rcr 

of aileron around c r 


tion 



natural angular frequency of wing in deflec- 
//Ch\ 

KVt) 

rrpb 2 

wing density parameter, where K = — ' 

M 

is the ratio of a mass of cylinder of air of 
a diameter equal to the chord of the wing to 
the mass of the wing, both taken for equal 
length along the span; this ratio may be 

expressed as K= 0 ,2i|.(b^/w) (p/p Q ) where W 

is the weight in pounds per foot span, b is 
in feet and p/p is the ratio of air density 
at altitude to that for normal standard air 

( _ M _ TT l\ 

K ~ ~ 5 V 



1+0 
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g a > Sr, Sh structural damping coqfficients (300 refer- 
ence 6 ) 


L 1> L 2' Ij 3> L V M l» M 2> M 3> quantities defined in table II 
and by equations (26) and (28) 


v/bco 


a 


QX = 
a 


flutter coefficient; that is, flutter speed 
divided by reference speed b<o a ' 



QpX — 
= 



v/hare .is the angular (flutter ) frequency and 



for case of bending- torsion (Hots that in the incom- 
pressible case (references ii and 6) u is replaced 
by 1/k.) 
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TABLE I.- VALUES OP f Q (M, w) = (t 0 ) R + 


lfo) R = | | J o(l) °° a 
(f o)i = “i J W J o{&) ain u du 






0.526 

1.053 

1.516 

1.754 

2.105 

2.506 

3.289 

4.211 

5.013 

6.266 

7.519 

9.569 

20.243 

40.486 

105.263 


M = i£ 


0.02107622 

.10786366 

.15423397 

.1759312$ 

.23564109 

.26927908 

.27020078 

.32453905 

.40860905 

.539882 

.64370304 


.75786234 

.938879 

.984312 

.997660 


-0.14998785 
- .21774161 
-.24646721 
-.27930193 
-.30623887 
-.29526037 

-!p3?7m 

-.47752247 

-.502746 

-.48895719 

-.144031451 

-.247350 

-.128388 

-.049910 



0.278 ■ 
.556 
1.111 
I.263 
1.684 
1.984 
2.525 

tag 

5.051 

6.313 

7.508 

9-921 

19.841 

39.683 

92.595 


-0.02589034 

.02529654 


.21539228 

.22643569 

.26380752 

.38184436 

.557683 

.678 777 

.79855505 


.844260 

.886947 

•933477 

.982907 

.995700 

.999200 


-0.08629977 

-.22399799 

-.33199853 

-.32383270 

-.36922184 

-.43392085 

-.51318210 

-.532604 

-.500827 

-.45177830 

-.387852 

-.338397 

-.266032 

-.138218 

-.069779 

-.029983 



20.00 

10.00 

5.00 
3.90 
3.10 
2.40 

2.00 
1.60 
l.l8 

:i 


0.196 

:?f 

1.00 

1.265 

1.634 

1.96H 

2.451 


3.323 

4.172 

5.028 

6.325 

m 

19.608 

39.216 

98.039 


01041793 - 
02790057 
13530140 
15895143 
.19688319 

59012790 

74934? 


8S1388 

923323 

945429 

967340 

991735 

997930 

999675 


M = f 


20.00 

10.00 

5.00 

3.10 

2.50 

1.90 

1.60 

1.30 


0.156 

.312 

.625 

1.008 

1.250 

1.645 

m 

3.324 

4.112 

5.040 

6.250 

7.U+9 

9.766 

19.531 

35.062 

78.125 


0.00827247 

-.03440806 

.07303312 

.15700194 

.59909796 

.71425408 

.839483 

.892157 

.926973 

.951528 

.965855 

.980047 

.994975 

.998738 

.999675 


0.05473581 

-.18976570 

-.33798972 

-.35747165 

-.44343786 


--■mm 

--■mr 

-.413781 

-.357217 

-.293258 

-.250038 

-.195428 

-.099425 

-.049930 

-.020000 


-0.07001922 

-.13439143 

-.52465775 

-.49133,028 

-.57452410 

-.60117841 

-.57737682 

-.52373582 

-.419702 

-.352880 

-.295089 

-.2421,10 

-.205300 

-.157913 

-.079738 

-.039962 

-.020000 
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TABLE I.- VALUES OP T Q (U, 3 ) = (f Q ) R + i(f 0 )j - Concluded 


37 

1 

]c 


Ml 

uj 

1 

E 

Mb 

f f o)i 

M = 2 

It = 5 

I 


-0.01548798 

-.00697070 

.00782279 

.21306136 

.41159770 

.60680594 

.72237641 

.79366285 

.885888 

.925636 

.946567 

.967381 

.975947 

.987392 

.996329 

.999333 

.99950 

-0.06537016 

-.08613197 

-.29541589 

-.59458735 

-.63550404 

-.59698731 

-.53566558 

-.47882893 

-.371721 

-.395330 

-.261128 

-.205798 

-.177347 

-.128996 

-.069843 

-.029983 

-.010000 

20.00 

10.00 

5.00 
4.20 
2.10 
1.70 
1.20 

1.00 

& 

•4° 

:1k. 

:?o° 

.06 

.02 

0.104 

.208 

& 

.992 

1.225 

1.736 

2.083 

2.480 

3.360 

4.167 

4.960 

6.127 

7.440 

10.417 

20.833 

34.722 

104.167 

- 0.1854996 

-.00207595 

-.15446620 

-.16645177 

.41122567 

.58072358 

.77380710 

.83908297 

.884607 

.936018 

.958080 

.970300 

.980474 

.986729 

.993215 

.998300 

.999383 

.999950 

-0.06011798 

-.12707440 

-.17734818 

-.35826002 

-.70315051 

-.65532307 

-.52771202 

-.45746891 

-.394511 

-.299626 

-.244532 

-.206750 

-.168271 

-.139932 

-.099645 

-.049960 

-.029983 

-.010000 

M = | 

M = i° 

3 

20.00 

10.00 

5.00 

4.80 

2.40 

1.90 

1.40 

1.20 

.96 

. 7.2 

$ 

.38 

.32 

.24 

.12 

.06 

.02 

0.119 

.238 

.476 

.496 

.992 

1.253 

1.701 

12ft 

3.307 

4.195 

4.960 

6.266 

7.440 

9.921 

19.841 

39.683 

119.048 

0.00671539 

.02216359 

-.05520172 

-.05938671 

.29706773 

.49182591 

.89153899 

.76519653 

.844299 

.909960 

.940834 

:?p6 

. 9 §l 697 

.989675 

.997403 

.999350 

.999950 

-0.04537548 
-.06784796 
-.25614005 
-.27691360 
-.65923330 
-.65290610 
-.57126881 
-.51684512 
-.436326 
-.341204 
-.28008 6 
-.234350 
-.187184 
-.158316 
-.119288 
-.059908 
-.029983 
-.010000 

20.00 

10.00 
5.00 

4.40 

2.20 

1.80 

1.30 

1.10 

.88 

.66 

'M 

.30 

.22 

.10 

.06 

.02 

0.110- 
.220 
. .440 
.500 
.999 
1.221 
I.69I 

l$i 

3.330 

4.227 

4.995 

6 . 10 S 

7*326 

9.990 

21.978 

56.630 

109.890 

0.00960890 

.05509080 

-.11086610 

-.12438944 

.37934955 

.801+20000 

.871331 

.926118 

:?! m 

•991595 

.998260 

.999367 

.999950 

-0.05304109 

-.08576219 

-.21422339 

-.31880518 

-.68956543 

-.65968944 

-.55534513 

-.49023486 

-.408825 

-.316656 

-.253427 

-.216005 
-.177806 
-.148727 
-.109495 
-.049950 
-.029983 ' 
-.010000 
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* V - 


TASLIt II.- VALUES OP FUNG TICK 3 USED II THE FLUTTER CALCULATIONS 


Tha axprs talons • cployoi In tb* ealonlatlon of this tsbls sr*t 



V ' h ♦ 4 i* + *l 


whars 



v = h - *i 

*Z’ ' H - *2 

V - - =z) - 4(4 * *i) 

81 ' » i['i Kf) “* E + » J °(j) “■ " * •’{!) * ln B ] 

% ,J (») “* 3 ' « J °£) * to * * '=© •“ 4 
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Figure 1.- Mach angle p . The disturbance at point (5,1]) moving forward with 
supersonic velocity v influences the angular region having half vertex angle 

sin-lf-. 



Figure 2.- Influence of impulse created at point (4,0) at time t - T on a point 
(x,y) fixed relative to (5,0) and moving with supersonic velocity v. (Observe 
that the disturbance influences the point (r j) only during the time interval 
2 " T l0 




Figs. 3,4 
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Figure 8.- Sketch showing that only disturbances created forward of the Mach 
angle region with vertex at £/, can affect (x,y). 


Tradliac 



Figure 4.- Sketch illustrating the three degrees of freedom h, a , and P of the 
oscillating airfoil. 
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Figure 5.- The flutter coefficient against center -of -gravity location for several positions of elastic axis 
and for three values of the frequency ratio. M - -^p-; h = 3.927. 
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Center-of- -gravity location, percent chord national advisory 

(a) Measured from leading edge. committee for aeronautics. 

(b) and (c) Measured from elastic axis. 

Figure 0.- The flutter coefficient- against center-of -gravity location for several positions of elastic axis 
and for three values of the frequency ratio. M = 7.854. 
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Figure 7.- The flutter coefficient against center-of -gravity location for several positions of elastic axis 
and for three values of the frequency ratio. M «= ; V-- 15.706. 


Fig. 7a-c 





Center -of -gravity location, percent chord 

(a) Measured from leading edge. advisory 

(b) and (c) Measured from elastic axis. COMMITTEE for AERONAUTICS. 


Figure 8.- The flutter coefficient against center -of -gravity location for several positions of elastic axis 
and for three values of the frequency ratio. M = p = 3.927 
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Center -of -gravity location, percent chord NATIONAL advisory 

(a) Measured from leading edge. COMMITTEE FO# AEBOHAUTICS. 

(bj and (c) Measured from elastic axis. 

Figure 9.- The flutter coefficient against center -of -gravity location for several positions of elastic axis 

and for three values of the frequency ratio. M = P = 7.854. 
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Center-of -gravity location, percent chord national advisory 

(a) Measured from leading edge. COMMITTEE FOR aeronautics. 

(b) and (c) Measured from elastic axis. 

Figure 10.- The flutter coefficient against center-of -gravity location for several positions of elastic axis 

and for three values of the frequency ratio. M = — ; R= 15.706. 
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(a) Measured from leading edge. 
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Figure 11.- The flutter coefficient against center -of-gravity location for several positions of elastic axis 

and for three values of the frequency ratio. M = 2; ^ = 3.927. 
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(a) Measured from leading edge. 
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figure 12.- The flutter coefficient against center -of -gravity location for several positions of elastic axis 

and for three values of the frequency ratio. M * 2; M- = 7.854. 




Fix. 12a-c NACA TN No. 1158 







Center-of -gravity location, percent chord 

(a) Measured from leading edge. 

(b) and (c) Measured from elastic axis. 
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Figure 13.- The flutter coefficient against center-of -gravity location for several positions of elastic axis 
and for three values of the frequency ratio. M = 2 ;^ = 15.708. 
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Center -of -gravity location, percent chord NATIONAL ADVISORY 

(a) Measured from leading edge. COMMITTEE F0I1 AERONAUTICS. 

(b) and (c) Measured from elastic axis. 

Figure 14.- The flutter coefficient against center-of -gravity location for several positions of elastic axis 
and for three values of the frequency ratio. M => 6; P - 3.927. 
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Center -of -gravity location, percent chord NATIONAL advisory 

(a) Measured from leading edge. COMMITTEE FOR aeronautics. 

(b) and (c) Measured from elastic axis. 

Figure 15.- The flutter coefficient against center-of -gravity location for several positions of elastic axis 

and for three values of the frequency ratio. M = 5; t 1 = 7.854. 
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Center -of -gravity location, percent chord 

(a) Measured from leading edge. 

(b) and (c) Measured from elastic axis. 
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Figure 16.- The flutter coefficient against c enter - of -gravity location for several positions of elastic axis 

a nd for three values of the frequency ratio. M - 5; n = 15.706. 
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Figure 17.- The flutter coefficient against Mach number for two locations of 

Wy. 

the center of gravity. Other parameters are = 0.707; a = 0; 

H = 7.864. a 


Fig. 18 
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figure 18.- Hie flutter coefficient against Mach number lor several values ol 

“h 

p, . Other parameters are ■; 0; x a » 0.2; * » 0. 
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Figure 19.- The flutter coefficient against frequency ratio for several values of 
M. Other parameters are a = 0; 2 ^ = 0.2; p. = 7.854. 




Fig. 20 
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Figure 20.- The flutter coefficient against frequency ratio for three values of 
x Q . Other parameters are M = a = 0; p = 7.854. 



Figure 21.“ Plot erf = 0. 
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(a) Flutter coefficient against axis- 05 ) Flutter coefficient against tor- (c) Flutter coefficient against wing- 

of -rotation position for several slonal damping coefficient for density parameter a for 

^ ^ = l5 ; 708 )- two values of x n (M =■ several values of Zq (M = 4£). 

Note that the range of x 0 y . . . , .. „ , . ■ 

narrows with increase in M ^ = 15.708). Negative damping The straight-line curve shown 

and disappears at M = 1.58 values are shown dashed and corresponds to M 3 = 0 

and Zq = 0.33. have no physical existence. = 0.327). 

Figure 22.- Curves for one -degree -of -freedom torsional instability. 
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